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The kinetics of a variety of transport-controlled processes can be reduced to the problem of de¬ 
termining the mean time needed to arrive at a given location for the first time, the so called mean 
first passage time (MFPT) problem. The occurrence of occasional large jumps or intermittent pat¬ 
terns combining various types of motion are known to outperform the standard random walk with 
respect to the MFPT, by reducing oversampling of space. Here we show that a regular but spatially 
heterogeneous random walk can significantly and universally enhance the search in any spatial di¬ 
mension. In a generic minimal model we consider a spherically symmetric system comprising two 
concentric regions with piece-wise constant diffusivity. The MFPT is analyzed under the constraint 
of conserved average dynamics, that is, the spatially averaged diffusivity is kept constant. Our 
analytical calculations and extensive numerical simulations demonstrate the existence of an optimal 
heterogeneity minimizing the MFPT to the target. We prove that the MFPT for a random walk 
is completely dominated by what we term direct trajectories towards the target and reveal a re¬ 
markable universality of the spatially heterogeneous search with respect to target size and system 
dimensionality. In contrast to intermittent strategies, which are most profitable in low spatial di¬ 
mensions, the spatially inhomogeneous search performs best in higher dimensions. Discussing our 
results alongside recent experiments on single particle tracking in living cells we argue that the 
observed spatial heterogeneity may be beneficial for cellular signaling processes. 

PACS numbers: 89.75.-k,82.70.Gg,83.10.Rs,05.40.-a 


I. INTRODUCTION 

_ Random search processes are ubiquitous in nature 
BE, ranging from the diffusive motion of regulatory 
molecules searching for their targets in living biological 
cells 0, Hi [mill i bacteria and animals searching for food 
by active motion 01, all the way to the spreading of 
epidemics and pandemics SI and computer algorithms 
in high-dimensional optimization problems [26] . The fact 
that the search strategy in these processes is to a large 
extent random reflects the incapability of the searcher 
to keep track of his past explorations at least over more 
than a certain period 0. During the years several dif¬ 
ferent search strategies have been studied in the litera¬ 
ture, including Brownian motion BE E E . spatio- 
temporally decoupled L evy flig hts (LFs) [20l fUl and cou¬ 
pled Levy walks (LWs) I35l I l6l| in which the searcher un¬ 
dergoes large re-allocations with a heavy-tailed length 
distribution either instantaneously (LFs) or with con¬ 
stant speed (LWs), as well as intermittent search pat¬ 
terns, in which the searcher combines different types 
of motion 0, for instance, three-dimensional and one¬ 
dimensional diffusion E - Ebl . EE , three-dimensional 
and two-dimensional diffusion ii , or diffusive and 
ballistic motion 0 0 HTTI [E. 

The efficiency of the search strategy is conventionally 
quantified via the mean first passage time (MFPT) de- 
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fined as the average time a random searcher needs to 
arrive at the target for the first time 0, E, US E, US ■ 
The physical principle underlying an improved search ef¬ 
ficiency is an optimized balance between the sampling 
of space on a scale much larger than the target and on 
the scale similar to or smaller than the target 0. More 
specifically, periods of less-compact exploration—for in¬ 
stance, diffusion in higher dimensions, Levy flights, or 
ballistic motion—aid towards bringing the searcher faster 
into the vicinity of the target. Concurrently a searcher 
in such a less compact search mode may thereby easily 
overshoot the target [32, Ej]. In contrast, compact explo¬ 
ration of space (for instance, diffusion in one dimension) 
is superior when it comes to hitting the target from close 
proximity but performs worse when it comes to the mo¬ 
tion on larger scales taking the searcher from its starting 
position into the target’s vicinity: typically frequent re¬ 
turns occur to the same location, a phenomenon referred 
to as oversampling. Mathematically this is connected to 
the recurrent nature of such compact random processes. 
The idea behind search optimization is to find an opti¬ 
mal balance of both more and less compact search modes 
in the given physical setting 0. For instance in the so- 
called facilitated diffusion model for the target search of 
regulatory proteins on DNA[E]23| the average duration 
of non-compact three-dimensional free diffusion is bal¬ 
anced with an optimal compact one-dimensional sliding 
regime along the DNA molecule. In intermittent search 
0ilE persistent ballistic excursions are balanced by 
compact Brownian phases. This optimization principle 
intuitively works better in lower dimensions, where a 
searcher performing a standard random walk oversam- 
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pies the space. Hence, the typically considered optimized 
strategies have the largest gain in low dimensions. 

In a variety of experimental situations the motion of 
a searcher is characterized by the same search strategy 
but is not translationally invariant. A typical example is 
a system in which the searcher performs a standard ran¬ 
dom walk but with a spatially varying rate of making its 
steps. This type of motion is actually abundant in biolog¬ 
ical cells, where experiments revealed a distinct spatial 
heterogeneity of the protein diffusivity jfiol-l&l]. Several 
aspects of such diffusion in heterogeneous media have al¬ 
ready been addressed [13, l63H65l | , but the generic FPT 
properties remain elusive, in particular, for quenched en¬ 
vironments. 

Here was ask the question whether spatial heterogene¬ 
ity is generic.ally detrimental for the efficiency of a ran¬ 
dom search process or whether it could even be beneficial. 
Could it even be true that proteins find their targets on 
the genome in the nucleus faster because their diffusivity 
landscape in the cell is heterogeneous? On the basis of 
exact results for the MFPT in one, two, and three di¬ 
mensions in a closed domain under various settings we 
here show that a spatially heterogeneous search can in¬ 
deed significantly enhance the rate of arrival at the tar¬ 
get. We explain the physical basis of this acceleration 
compared to a homogeneous search process and quan¬ 
tify an optimal heterogeneity , which minimizes the MFPT 
to the target. Furthermore we show that heterogeneity 
can be generically beneficial in a random system and is 
thus a robust means of enhancing the search kinetics. 
The optimal heterogeneous search rests on the remark¬ 
able observation that the MFPT is completely dominated 
by those trajectories heading directly towards the target. 
We prove that the MFPT for the heterogeneous system 
can be exactly described with the results of a standard 
random walk. We compare our theoretical findings to 
recent experiments on single particle tracking in living 
cells, which are indeed in line with the requirements for 
enhanced search. 

The paper is organized as follows. Section II introduces 
our minimal model for heterogeneous search processes. In 
Sec. Ill we briefly summarize our main general results, 
which hold irrespective of the dimension (d = 1, 2, or 3). 
Section IV is devoted to the analysis of the most general 
situation with a specific starting point and position of the 
interface. In Sec. V we focus on the Global MFPT, that 
is, the MFPT averaged over the initial position. In Sec¬ 
tion VI we analyze a system with a random position of 
the interface and optimize the MFPT averaged over the 
interface position. In Section VII we address the Global 
MFPT in systems with a random position of the inter¬ 
face. Throughout we discuss our results in a biophysical 
context motivated by recent experimental findings. Fi¬ 
nally, we conclude by discussing the implications of our 
results for more general spatially heterogeneous systems. 


II. MINIMAL MODEL FOR SPATIALLY 
HETEROGENEOUS RANDOM SEARCH 

We focus on the simplest scenario of a spatially hetero¬ 
geneous system. Even for this minimal model the analy¬ 
sis turns out to be challenging and our exact results reveal 
a rich behavior with several a priori surprising features. 
We consider a spherically symmetric system in dimen¬ 
sions d = 1,2, and 3 with a perfectly absorbing target of 
radius a located in the center (Fig.[T|). The outer bound¬ 
ary at radius R is taken to be perfectly reflecting. The 
system consists of two domains with uniform diffusivi- 
ties denoted by D\ and D 2 in the interior and exterior 
domains, respectively. The interface between these do¬ 
mains is located at r\. The microscopic picture we are 
considering corresponds to the ’kinetic’ interpretation of 
the Langevin or corresponding Fokker-Planck equations. 
In particular we assume that the dynamics obey the 
fluctuation-dissipation relation and in the steady state 
agree with the results of equilibrium statistical mechan¬ 
ics iod • 

In the biological context we consider that the system 
is in contact with a heath bath at constant and uni¬ 
form temperature T. The signaling proteins diffuse in 
a medium comprising water and numerous other par¬ 
ticles, such as other biomacromolecules or cellular or¬ 
ganelles. The remaining particles, which we briefly call 
crowders, are not uniformly distributed across the cell- 
their identity and relative concentrations differs within 
the nucleus and the cytoplasm and can also show vari¬ 
ations across the cytoplasm. The proteins hence expe¬ 
rience a spatially varying friction T(r), which originates 
from spatial variations in the long-range hydrodynamic 
coupling to the motion of the crowders, which is in turn 
mediated by the solvent [(37| [68|. The proteins thus move 
under the influence of a position dependent diffusion coef¬ 
ficient D(r) = 2fcsTr(r) and a fluctuation-induced ther¬ 
mal drift F(r) ~ fcBTVr(r) (see j69[ for details). The 
vital role of such hydrodynamic interactions in the cell 
cytoplasm was demonstrated in [70]. 

Because of the spherical symmetry of the problem we 
can reduce the analysis to the radial coordinate alone. 
That is, we only trace the projection of the motion of the 
searcher onto the radial coordinate) and therefore start 
with a discrete space-time nearest neighbor random walk 
in-between thin concentric spherical shells of equal width 
A R = Ri + i — Ri as depicted in Fig. [T] The shell i denotes 
the region between surfaces with radii Ri -1 and Ri. We 
assume that the hopping rates between shells n(z —> j) 
obey detailed balance 

p(i)Il(i —* + 1) = p{i + l)n(* + 1 — > i). (la) 

Here p(i) denotes the probability distribution of finding 
the particle in shell i. The hopping rates n are given as 
the product of the intrinsic rate 2D(i)/AR 2 and q[i ), the 
probability to jump from shell i to shell i +1 (and 1 — q(i) 
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FIG. 1: Schematic of the model system: a) A spherical target with radius a is placed in the center of a spherical domain 
of radius R. The free space between the radii a and R is divided into two regions denoted by subscripts. The inner region 
is bounded by a shell at radius n. The outer region ranges from r\ to the reflective boundary at R. Initially, the particle’s 
starting position is uniformly distributed over the surface of the sphere with radius r o. b) Microscopic picture of the problem 
starting from a discrete random walk between spherical shells. The hopping rates are assumed to obey detailed balance and 
the interface position is chosen to be placed symmetrically between two concentric spherical surfaces. 


for jumps in the other direction), 

n(i^i + l) = ^q(i). (lb) 

Here D(i ) is the arithmetic mean of the diffusivity in 
shell i. The rate q(i) can be derived as follows. A ran¬ 
dom walker located in shell i at time t can either move 
to shell i + 1 with probability q(i) or to shell i — 1 with 
probability 1 — q(i). Due to the isotropic motion of the 
random walker these probabilities are proportional to the 
respective surface areas of the bounding d-dimensional 
spherical surfaces at Ri-i and It ,. That is, q{i) ~ R d ^ 1 


and 1 — q{i) ~ R-fli ■ The proportionality constant is 
readily obtained from the normalization condition lead¬ 
ing to 


q{i) = 


R" 


d -1 


Rt l + RtV 


(lc) 


d ^ Therefore, 


and thus 1 - q{i) = R d Zi/{R d ~ + R d Zi)- 
while the random walker moves in all directions (radial, 
azimuthal, or polar) we can project its motion on the ra¬ 
dial coordinate only. We assume that the interface is lo¬ 
cated between two concentric shells leading to a continu- 
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ous steady state probability density profile. The searcher 
starts at t = 0 uniformly distributed over the surface of 
a d-sphere with radius r o, as sketched in Fig. [T^,) . 

In our analytical calculations we model the system in 
terms of the probability density function p(r, f |ro) to find 
the particle at radius r at time t after starting from radius 
ro at t = 0. p(r,t\ro) obeys the radial diffusion equation 


dp(r,t\r 0 ) 

dt 


1 d 

r d ~l Q r 



P(r,t\r 0 ) 


(2a) 


with piece-wise constant diffusivity 


D(t) 


D i, a < r < r\ 
D 2 , n <r < R 


(2b) 


We assume that the target surface at radius a is perfectly 
absorbing, 


p(a,t\r o ) = 0, (2c) 


to determine the first passage behavior. At the outer 
radius R we use the reflecting boundary condition 


dp{r,t\r 0 ) 

dr 


= 0 . 

r—R 


(2d) 


These boundary conditions are complemented with join¬ 
ing conditions at n by requiring the continuity of the 
probability density and the flux, which follow from our 
microscopic picture. 

To quantify our model system we introduce the ratio 

Di , , 

V = 77 - (3a) 


of the inner and outer diffusivities. Moreover, we demand 
that the spatially averaged diffusivity 


D = 


R d - a d 


r d ~ l D(r)dr 


(3b) 


remains constant for varying D\ and Z? 2 - Without such 
a constraint the problem of finding an optimal p, which 
minimizes the MFPT is ill-posed and has a trivial so¬ 
lution ip = 00 . More importantly, we want to compare 
the search efficiency as a function of the degree of het¬ 
erogeneity, where the overall intensity of the dynamics is 
conserved. Returning to our microscopic picture of a sig¬ 
naling protein searching for its target in the nucleus, the 
heterogeneous diffusivity is due to spatial variations in 
the long-range hydrodynamic coupling to the motion of 
the crowders. Their identity and relative concentration 
in the cell varies in space, but the effect is mediated by 
thermal fluctuations in the solvent at a constant temper¬ 
ature. The constraint in Eq. (13bl) then corresponds to a 
redistribution of the crowders at constant temperature, 
cell volume and numbers of the various crowders, which 
would not affect the spatially averaged diffusivity. 


Under the constraint © of constant spatially aver¬ 
aged diffusivity we obtain for any given p and n that 


D 1 


pD 

(<p - i)x(n) +1 ’ 


d 2 


D 

(<p - i)x(n) +1 ’ 


(4a) 


where we introduced the hypervolume ratio 


x(n) 



R d -a d ' 


(4b) 


of the inner versus the entire domain excluding the target 
volume. To solve Eq. (I2al) we take a Laplace transform 
in time and the obtained Bessel-type equation is solved 
exactly as shown in Sec. IV. The MFPT T(ro) for the 
particle to reach the target surface at r = a is obtained 
from the Laplace transformed flux into the target 


j(r 0 ,s) = fldDxa 1 


d -1 dP(r,r 0 ,s ) 


dr 


(5a) 


via the relation 


TM 


djai.ro, s) 


ds 


s—0 


(5b) 


The angular prefactor fid = 1 for d = 1, fid = for 
d = 2 and fid = 47t for d = 3. We treat the degree 
of heterogeneity p as an adjustable parameter at a fixed 
value of ri and seek for an optimal value minimizing the 
MFPT. The optimal heterogeneity, which we denote with 
an asterisk, is thus obtained by extremizing T a (ro) with 
respect to p. 


III. SUMMARY OF THE MAIN RESULTS 


Since R in combination with the diffusivity D sets the 
absolute time scale we can express, without loss of gener¬ 
ality, time in units of R 2 /D, set D = 1 and focus on the 
problem in a unit sphere. We introduce dimensionless 
spatial units x a = a/R, X\ = r\/R and Xo = xq/R. For 
the sake of completeness, we retain the explicit R and D 
dependence in the prefactors. 

Our first main result represents the fact that the 
MFPT to the target in the inhomogeneous system in di¬ 
mension d = 1, 2 and 3, can be expressed exactly in 
terms of the corresponding MFPT T°(a;o;IIi) in a ho¬ 
mogeneous system with diffusivity Di , with i = 1 or 2. 
Remarkably the MFPT is thus exactly equal to 


T0r o ) 


T 0 (x 0 ;Hi), x 0 < x\ 

T 0 ixi-D 1 ) + T° xl ixo;D 2 ), x 0 > x\ 


In the second line the argument xi stands for the re¬ 
lease of the particle at the interface, and the index xi 
of the last term is used to indicate that this term mea¬ 
sures the first passage to the interface at x\. That is, 
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FIG. 2: Schematic of the equivalence of MFPTs in inhomogeneous and homogeneous systems in the case of a) a searcher 
starting in the inner region and b) a searcher starting in the outer region. The radius x\ of the initial particle position is shown 
by the thin dashed circle. We show direct trajectories as full lines. Each panel also contains an indirect trajectory (dashed 
line) which leads the particle into the outer region of the system. 


in this case when the particle starts in the inner region 
with diffusivity D\ Eq. © reveals that the MFPT of the 
heterogeneous system is equal to that of a homogeneous 
system with diffusivity D\ everywhere and is indepen¬ 
dent of the position of the interface. Conversely, if the 
searcher starts in the exterior region with diffusivity D 2 
the MFPT contains two contributions: (i) the MFPT 
from tq to ri in a homogeneous system with diffusivity 
D 2 and (ii) the MFPT from r\ to a in a homogeneous sys¬ 
tem with diffusivity Di, as shown schematically in Fig. [2] 
Eq. © is exact and independent of the choice for D\p 
and thus holds for an arbitrary set of diffusivities and 
even if D\ = D 2 . In other words, it is not a consequence 
of a conserved D. 

The additivity principle of the individual MFPTs in 
Eq. © is only possible if the excursions of the searcher 
in the directions away from the target are statistically 
insignificant. We would expect that some trajectories 
starting in the inner region will carry the searcher into 
the outer region with diffusivity £>2 before they eventu¬ 
ally cross the interface and reach the target by moving 
through the inner region with diffusivity D\. Such tra¬ 
jectories will obviously be different in the heterogeneous 
system in comparison to a homogeneous system with dif¬ 
fusivity D\ everywhere, This appears to contradict the 
complete independence of D 2 of the MFPT in Eq. © 
for trajectories with xo < x\. This observation can be 
explained by the dominance of direct trajectories , whose 
occupation fraction outside the starting radius is statis¬ 
tically insignificant: The MFPT for a standard random 
walk in dimensions d = 1,2, and 3 is completely domi¬ 
nated by direct trajectories. As such, the MFPT is really 
a measure for the efficiency of the fast trajectories. 

Our second main result demonstrates that for xo > x\ 
a finite optimal heterogeneity exists at given interface 


position and is given by 


ip*(x 0 ) 


1 - X(si) x TQ(x i; 1) 
X(ah) T° t (a; 0 ; 1)' 


(7) 


For this value the MFPT T(a;o) attains a minimum. 
Hence, the optimal heterogeneity is completely deter¬ 
mined by the volume fractions and the MFPT properties 
and hence strictly by the direct trajectories. As above, 
the index xj in the MFPT T°(a;i indicates the first pas¬ 
sage to the interface, while without this index the MFPT 
T° quantifies the first passage to the target at x a ■ The 
explicit results for (p*(x 0 ) are shown in Fig. 0] and are 
discussed in detail in Sec. IV. 

Often one is interested in the MFPT averaged over an 
ensemble of starting positions, the Global MFPT T. As 
before, it can be shown that an optimal heterogeneity 
exists for any interface position and is universally given 
by 


T 


f i - xQi) 
V x{xi) 


T°(xi;l )-fZ 

Zo^T°(zo;l)' 

dx 0 


Zo^TO^o;!)' 

dx 0 


1/2 

• ( 8 ) 


Similar to the general case Tp* is again proportional to 
(x(xi) -1 — l) 1 / 2 but here the corresponding MFPTs in 
the second factor are reduced by a spatially averaged 
change of the MFPT with respect to the starting posi¬ 
tion. The optimal heterogeneity for the Global MFPT is 
shown in Fig. 0 and discussed in detail in Sec. V. 

In a setting when the interface position is random and 
uniformly distributed we are interested in the MFPT 
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from a given starting position averaged over the inter¬ 
face position. A measurable quantity for this scenario for 
an ensemble of random-interface systems is the MFPT 
{T a (x 0 )}, where the curly brackets denote an average 


over the interface positions x\. Explicit results for di¬ 
mensions d = 1,2, and 3 are given in Sec. VI. Solving for 
the optimal heterogeneity we obtain 


M* = 


t ° m -K )/ 


X ( 1 — 


d + 1 


^-T°(x 0 ) 

dx i 


1/2 


dx i 


d+1 rpO 


V 


T> 0 )- 1 + 


1 


d + 1 


- x. 




dx i 


(9) 


The optimal heterogeneity for the MFPT averaged over the random interface position is shown in Fig. [Gj 

Finally, we compute the Global MFPT averaged over the position of the interface, {T}, whose explicit results are 
given in Sec. VI. Also here an optimal strategy can be identified as 


(¥ 


m* = 


T° - 


/ 1 / . j 1 —{— d Xj 1 d __n . , , 

yji-xOx, i 


1/2 


1 T _ 


xi 4 - [1 +M 4-x°..fe 


l-®o 

I- 


dx i 


L Xl Uj ' X 'Xl 


( 10 ) 


The optimal heterogeneity for the Global MFPT av¬ 
eraged over the random interface position is shown in 
Fig. 0 

Eq. (0 represents a rigorous proof that direct trajec¬ 
tories dominate the MFPT for Brownian motion. In ad¬ 
dition, the heterogeneity does not affect the fraction of 
direct versus indirect trajectories but only their respec¬ 
tive durations. Due to the fact, that indirect trajectories 
are statistically insignificant, we can in principle make 
them arbitrarily slow if we start in the inner region. But 
as we let the inner diffusivity go to infinity (and hence the 
outer one to zero) we are simultaneously slowing down 
the arrivals to the interface if starting from the outer 
region. The physical principle underlying the accelera¬ 
tion of search kinetics is: The optimal heterogeneity cor¬ 


responds to an improved balance between the MFPT to 
reach the interface and the one to reach the target from 
the interface. 


IV. MEAN FIRST PASSAGE TIME FOR FIXED 
INITIAL AND INTERFACE POSITIONS 

Here we present the mathematical derivation and the 
explicit results for the situation with a specific initial con¬ 
dition ro and interface position n. Eq. Eal) is solved by 
Laplace transformation, and the resulting Green’s func¬ 
tion with the boundary conditions (Eel) and (I2dl) reads 


P(r, s|r 0 ) = 


(rr 0 y 


C v (S ir ,S ia ) 


T> v (S\ro : Sdn) , 1 C v (S\ro, Sin) 


& d Di T> v {Sia,Sxn) 1 C^S^Sxn) 


C v . 1 (S 2 r 1 ,S 2 R) JipV v (S 2 ri,S 2 Ry 
■D„(S 2 r 0 ,S 2 R) 


C v - l {S 2 n,S 2 R) y/tpV v (S 2 n,S 2 R) l nSiT>u{S 2 ri, S 2 R)C„-i(S 2 ri, S 2 R) ’ 


a < ro < n 

n <r 0 < R 


I1M 


where we introduced the abbreviation S\ t2 = \JsfD\ i2 
and the auxiliary functions 


V u (z\,z 2 ) = I v (z\)K v - 1 {z 2 ) + K v {z\)I v -i(z 2 ) 
C u (zi,z 2 ) = I i/ (z 1 )K v (z 2 ) - I„(z 2 )K v (z 1 ). (lib) 


Here the I v {z) and K„(z) denote the modified Bessel 
functions of order v = 1 — d/2 of the first and second kind, 
respectively. The Laplace transformed first passage time 
density is obtained from the flux (l5al) into the target, 
and the MFPT then follows from relation (I5bl) . In case 
of d = 1 the target size only enters the problem by deter- 
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mining the width of the interval. Using T> v (z , z) = \j z it 
can be shown that Eq. (1 11 all reduces to the ordinary ex¬ 
pression for regular diffusion given in Ref. sa for r\ = R 
and p = 1. 

The MFPT can be written exactly in terms of the ex¬ 
pressions for a homogeneous diffusion process in Eq. © 
(compare Refs. 0,01), and we obtain 


iT°(x 0 ;D) 
2 T°(x 0 ;D) 
3 T°(x 0 ;D) 


_R2 

W' 

4 D 
R 2 

W 


[2(x 0 - x a ) + x 2 - a$, 
2 log 


( x o 

U 


, 2 2 
X n Xr\ 


Xq X a 2 2 


(12a) 

(12b) 

(12c) 


Here the left index denotes the dimensionality of the sys¬ 
tem. Our results show excellent agreement with numeri¬ 
cal simulations, as demonstrated in the Appendix. Plug¬ 
ging the above expressions (1121) into Eq. (O we can com¬ 
pare the search efficiency with respect to a homogeneous 
random walk by introducing the dimensionless ratio 


d(x 0 ) 


T(® 0 ) 

T°(x 0 ) ’ 


(13) 


The corresponding results for dimensions d = 1, 2, and 3 
are shown in Fig. [3] 

The qualitative behavior of the MFPT with respect 
to p —that is, the degree of the heterogeneity—depends 
on the starting position relative to the interface. If 
the initial position lies in the inner region 9(xg ) decays 
monotonically and saturates at a finite asymptotic value, 
9(x o) — > x( x i)- Hence, in this case an optimal strat¬ 
egy does not exists and the best search performance is 
achieved for large diffusivities in the inner region. The 
lower bound on 9(x o) is set by the volume fraction of 
the inner region, which sets a bound on the ratio Di/D. 
This result is yet another manifestation of the fact that 
the MFPT is completely dominated by direct trajecto¬ 
ries. 

Conversely, if the searcher starts in the outer region 
an optimal strategy exists. To understand the existence 
and meaning of the optimal heterogeneity we perform 
a power series expansion of 9{x o). We find the scaling 
9{x o) ~ 1/p as p — >0, which is due to the fact that in 
this regime D± ~ p 7 which dominates the MFPT in this 
regime. In the other limit as p — > oc we find 9(xo) ~ p 
because here Di ~ p and the rate determining step is the 
arrival at the interface. We can understand the optimal 
heterogeneity as a beneficial balance between the rate of 
arriving at the interface from the starting position and 
the rate to find the target if starting from the interface. 


In contrast, too large values of p prolong the time 
to reach the interface and cannot be compensated by a 
faster arrival from the interface towards the target. The 
optimal heterogeneity as a function of the starting posi¬ 
tion is shown in Fig. |4ji)-c) and reveals the divergence as 
xg —> x\: this point corresponds to the disappearance of 
an optimal strategy. As xg gradually approaches unity, 
p* continuously decreases towards a plateau, which sug¬ 
gests that to reach the target from the interface becomes 
the rate determining step. Moreover, p* decreases with 
increasing target size, because the inner region becomes 
smaller and thus allows a larger Z ?2 hi the optimal sce¬ 
nario. The overall gain of an optimal search with re¬ 
spect to a standard random walk is shown in Fig. 0]I)- 
f). Only in the domain xg > x\ an optimal heterogene¬ 
ity exists, therefore we omitted the region Xg < x\. As 
mentioned before, we observe the monotone convergence 
9{x o) —> x( x i) from above as p —> oo and hence the 
heterogeneous search always outperforms the standard 
Brownian search for p > 1. The highest gain is therefore 
set by the volume fraction of the inner region which be¬ 
comes arbitrarily small as x\ — > x a . For Xg > x\ we find 
that the gain is largest for starting positions near the in¬ 
terface and when the interface is not too close to the outer 
boundary, where the system essentially approaches the 
homogeneous limit. The variations are larger in higher 
dimensions, which is, of course, connected with the pro¬ 
nounced spatial oversampling in lower d. 

Different from conventional strategies—intermittent or 
Levy-stable processes, which affect the compactness of 
exploring space—heterogeneous search is more profitable 
in higher dimensions (see Fig. 0Ji)-f)). Because hetero¬ 
geneity acts by enhancing/retarding the local dynamics 
and does not affect spatial oversampling it is intuitive 
that it performs better for non-compact exploration of 
space. Both the existence as well as the gain of an opti¬ 
mally heterogeneous search are thus a direct consequence 
of direct trajectories dominating the MFPT. 


V. GLOBAL MEAN FIRST PASSAGE TIME 
FOR FIXED INTERFACE POSITION 


The Global MFPT is obtained by direct averaging of 
Eq. ([6]) over the initial position xg , 

d r 1 

T = _ d / XQ~ 1 T(xg)dxg. (14) 

^ X a Jx a 

The exact expressions for the Global MFPT in the vari¬ 
ous dimensions read 












FIG. 3: Ratio 9(x o) = T(*o)/T°(a;o) as function of p = D 1 /D 2 in dimension a) d = 1, b) d = 2, and c) d = 3 for initial radius 
xo = 0.3 (dashed lines) and xo = 0.8 (full lines), respectively. In b) and c) we plotted the results for a target size x a = 0.1. 
The full black lines denote the ~ 1 /tp and ~ tp scaling, respectively, and the dashed black line corresponds to ~ x{ x i)- 
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FIG. 4: Optimal heterogeneity as function of the starting position xo for three interface positions xi = 0.25 (blue), xi = 0.5 
(red), and xi = 0.75 (orange) in dimension a) d = 1, b) d = 2, and c) d = 3. The target sizes in b) and c) are x a =0.1 (full 
lines) and x a =0.2 (dashed lines). Ratio 9{x 0 ) = T(a;o)/T°(xo) for the optimal heterogeneity tp* in dimension d) d = 1, e) 
d = 2, and f) d = 3. In e) and f) the target size is x a = 0.1. 
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These are to be compared with their homogeneous coun¬ 
terparts 


T o R 2 

lT = W 
2 1 = 


A 2 


4Zl(l-x 2 ) 2 


bT = -=- 


+ 2 X 2 - y - 2l0g(x o ) 

R 2 


3H(l-x 3 ) 2 
x 


9x„ 3 

5 + “ 5 


(16a) 


(16b) 


(16c) 


in d = 1, 2, and 3, respectively. As before we introduce 
the dimensionless enhancement ratio 


T°' 


(17) 


The result are shown in Fig. [5] In this case we are ef¬ 
fectively considering a weighted average of the results 
presented in the previous Sec. IV. Noticing that most 
of the volume of a d-sphcre is increasingly concentrated 


that the results will be more prominently influenced by 
the features of trajectories starting further away from the 
target. An interesting question will therefore be whether 
an optimal heterogeneity exists for all interface positions 
and dimensions. An expansion of T in a power series of 
ip reveals that 9 ~ 1/tp as tp — >0 and 9 ~ tp as p —> oo 
for all positions of the interface (see Fig. [5]a)-c)). Thus, 
an optimal heterogeneity always exists for all d. From 
Fig. [5] a)-c) we observed that the dependence of the rel¬ 
ative gain compared to the 6 on the interface position in 
the limit —>■ 0 is very weak in all dimensions as long 
as it is not too close to the external boundary. In the 
case of large p somewhat larger variations are observed 
for dimensions d = 1 and 2. In d = 3 the dependence 
on the interface position is largest near the optimal value 
for p but very weak elsewhere. Hence, a control of the 
target location dynamics by adjusting x\ is only efficient 
near the optimal point Tp = Tp*. This is observed from 
the respective scaling 0 ~ 1/p and 9 ~ ip as p — > 0 and 
p —> oo as explained in Sec. IV. Overall the gain with 
respect to the homogeneous random walk is larger for 
higher dimensions, which has the same origin as in the 
general case discussed above, however here the additional 
effect of averaging over the initial position enters. The 
optimal heterogeneity has the exact results 


IP = 


2P = 


3P = 


dimensions we 

anticipate 

1 

V / 3( 1 — xi) + x 2 


1 — xi 
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4 + Xi 2 + X 2 + 2 2 log f ) 
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1/2 


—41og(xi) — 3 + 4x 2 — Xj 
1 - x(®l) + *«) - (! + *l)] + X ^ x l ~ x a)' 


(18a) 

(18b) 


xi x i) 


(1 — xi) 3 (5 + xi [6 + xi(3 + xi)]) 


1/2 


(18c) 
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The behavior of Eqs. m is shown in Fig. [5j;)-f). In 
higher dimensions the optimal heterogeneity shows a 
non-monotonic behavior with respect to the interface po¬ 
sition and increases upon approaching the target or the 
outer boundary. Simultaneously, the overall dependence 
of the Global MFPT on ip vanishes in these limits—see 
the inset of Fig. |SJ)e)-f)— corresponding to the situation 
when the system is no longer heterogeneous. These re¬ 
sults can be rationalized by the fact that in the limit 
xi —> x a the ratio x{ x i) becomes negligible and hence a 
higher Di is allowed without slowing down the dynamics 
of reaching the interface from the external region. Con¬ 
versely, as x\ —> 1 also x( x i) V 1 and smaller D 2 are 
allowed because the rate limiting step is hitting the tar¬ 
get from the inner region. In both limits, however, the 
overall enhancement effect with respect to a standard 
random walk becomes negligible. Away from these limits 
the gain of the optimal heterogeneity is larger for higher 
dimensions and can be remarkably large for intermedi¬ 
ate interface positions (see the inset of Fig. [5])e)-f) and 
increases with decreasing target size. 

In a biological context, the present setting is relevant 
for signaling proteins searching for their target in the nu¬ 
cleus when the proteins are initially uniformly distributed 
throughout the cell cytoplasm. Recent experiments re¬ 


vealed a significant heterogeneity in the spatial depen¬ 
dence of the protein diffusion coefficient across the cell, 
with a faster diffusivity near the nucleus [(3(1 L61| . Such a 
spatial heterogeneity could therefore be beneficial for the 
cell by accelerating the dynamics of signaling molecules. 


VI. MEAN FIRST PASSAGE TIME IN A 
RANDOM HETEROGENEOUS SYSTEM FOR 
FIXED INITIAL POSITION 

We now address the MFPT problem when the interface 
position is random in a given realization and uniformly 
distributed over the radial domain. Specifically, we are 
interested in the MFPT of a particle starting at Xq aver¬ 
aged over the interface position xi, 

1 f 1 

{T}(x 0 ) = —— / T(x 0 )dxL (19) 

-L X a J Xa 

Experimentally, this would correspond to measuring an 
ensemble of systems with random value of xi. For dimen¬ 
sions d = 1, 2, and 3 the MFPT {T}(;ro) has the explicit 
form 


i{T}(*o) 

2 {T}(zo) 
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(20b) 


(20c) 


The gain compared to the homogeneous random walk, 


{6»(x 0 )} 


{T}(x 0 ) 

T(z 0 ) ’ 


( 21 ) 


is shown in Fig.(nji)-c) and reveals a somewhat less sharp 
optimal heterogeneity as compared to the Global MFPT. 


The limiting behavior for small and large ip is the same as 
for the previous cases but here there is a wider range of 
ip values producing a comparable gain. We also observe 
a stronger dependence on the starting position as ip —> 
00 and a very weak dependence for —>■ 0, which is 
yet another consequence of direct trajectories dominating 
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FIG. 5: Ratio 6 = T/T° as function of ip = D 1 /D 2 for different dimensions and interface positions xi = 0.4 (blue) and xi = 0.6 
(orange). The target size in panels b) and c) is x a = 0.1 (full lines) and x a = 0.2 (dashed lines). d)-f) Optimal heterogeneity 
for different dimensions as a function of the interface position xi. In panels e) and f) the target size is x a =0.05 (blue), x a =0.1 
(red), and x a =0.2 (orange). Insets in panels d)-f): 8 for the optimal heterogeneity Tp*. 


the MFPT. The optimal heterogeneity is obtained from Eq. m and reads 

I 


iM* 

2M* 

aM* 


J_ f 3(2 — z 0 )(2 — Xq[ 2 — Xp]) \ 1/2 
xq V 4 - 3x 0 / 

( 21og(x 0 /x a ) - 3[Qi(x 0 ) - Q 1 (x a )] \ 1/2 
V 2xl\og(x 0 /x a ) + Q 2 (xo,Xa) ) 

( 24log(xp/x a ) + Q 3 (x 0 ) - Q 3 (zaA 1/2 
\ 24x3 log(x 0 /x Q ) - < 34 (^ 0 , Xa) ) 


Here we introduced the auxiliary functions 


Qi(y) 

Q2(y,z) 

Qs(y) 

Q±(y,z) 


y 3 + y 


r 

35 




18 


+ y 2 (9-y[10-y 3 }), 


y 3 (2[l + y 2 z\ — y 3 ) + z 3 (2 + y 3 ) 





(22a) 

(22b) 

(22c) 

(23a) 

(23b) 

(23c) 

(23d) 




































12 


The results for various dimensions are depicted in 
Fig. EJi)®). Accordingly an optimal heterogeneity al¬ 
ways exists. Note that the dependence of {</?}* on the 
target size in d = 2 and 3 becomes reversed upon increas¬ 
ing the initial separation to the target. When starting 
near the target the initial position will on average lie in 
the inner region in an ensemble of realizations of the in¬ 
terface position. The optimal heterogeneity will thus on 
average correspond to a very fast diffusion in the inner 
region. For small xq the probability of starting in the 
inner region, (1 — x 0 )/(l — x a ), will be lower for smaller 
targets and hence {<^}* will be larger accordingly. Con¬ 
versely, if starting further away from the target xq will on 
average lie in the outer region and the search time will 
be more strongly influenced by the rate of arriving at 
the interface in each realization. An optimal heterogene¬ 
ity will therefore correspond to a smaller asymmetry of 
diffusivities in the inner and outer regions. More specifi¬ 
cally, since the probability of starting in the outer region, 
(a;o—x Q )/(l—x a ), will be lower for larger targets the value 
of {(/?}* will accordingly be smaller. The gain of the op¬ 
timal heterogeneity is shown in the insets of Fig. [HJl)-®) 
and reveals a significant improvement with respect to the 
standard random walk in every dimension. The depen¬ 
dence on target size has again a non-monotonic behav¬ 
ior, which follows from the argument presented above. 
Hence, even in a system with a random and quenched 
position of the interface a spatially heterogeneous search 


process will on average strongly outperform the standard 
random walk in every dimension. 

In a biological context the present setting is important 
for the experimentally observed quenched spatial disor¬ 
der revealed by single-particle tracking [6ll . 16^ 1651 . p7ll473| . 
More generally, the ideas can also be extended to the dy¬ 
namics in complex disordered systems Bill and Sinai- 
type diffusion [76|. Our results show that quenched spa¬ 
tial heterogeneity can robustly enhance random search 
processes even in a disordered setting. This robustness 
could eventually be exploited in search strategies, be¬ 
cause it requires less prior knowledge about the location 
of the target. 


VII. GLOBAL MEAN FIRST PASSAGE TIME 
IN A RANDOM SYSTEM 

This section completes our study, addressing the 
Global MFPT in a random system, 

{T} = -±— / Tdxi. (24) 

1 X a Jx a 

As far as blind spatially heterogeneous search is con¬ 
cerned this is the most robust setting. It can be shown 
that the exact results for {T} have the form 
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(25b) 
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The relative gain 


m = 


{T} 


(26) 


with respect to the Global MFPT for the standard random walk is shown in Fig. [7] The limiting scaling as tp —> 0 
and tp —y oo remains unchanged and the dependence of the overall gain on the target size becomes significant for large 
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FIG. 6: Ratio {#} = {T}/T(xo) as function of ip = D 1 /D 2 for different dimensions and *0 = 0.4 (blue) and xo = 0.6 (orange). 
The target size in panels b) and c) is x a = 0.1 (full lines) and x a = 0.2 (dashed lines). d)-f) Optimal heterogeneity for different 
dimensions as a function of the starting position xo■ In panels e) and f) the target size is x a = 0.05 (blue), x a = 0.1 (red), and 
x a = 0.2 (orange). Inset in panels d)-f): {0} for the optimal heterogeneity {y>*}. 


ip. The optimal heterogeneity in the different spatial dimensions is obtained in the form 


iM* 

aM* 

3M* 


= A, 


— / —127 - 84log(xa) + x a (126 + x a [84 - x a {98 + 3x a (7 - 2x a [7 - x£])}]) 
° v 16 - 420x3 log(x a ) - x 2(336 - x a [245 + 3x2{28 - 3x2}]) 

135/x a - 148 + x 2(135 - x a [135 - x^{9 - x^}]) + 180 log(x 0 )\ 1/2 
5 + x3(210-x Q [243-x2{30-2x3}]) + 180x3 1og(x Q ) J ' 


1/2 


(27a) 

(27b) 

(27c) 


The results are plotted in Fig.[7ji). It can be shown that 
the results in the d = 2 and d = 3 cases converge to the 
the result obtained in d = 1 in the limit x a —> 1. This 
is expected since the system effectively becomes one di¬ 
mensional when the ratio annulus thickness-to-curvature 
approaches zero. In the other limit x a —> 0 we find the 
diverging optimal heterogeneity (see Fig. 0))): 


in d = 3. The gain of the optimal heterogeneity is shown 
in Fig. [7}:). In d = 1 the optimal gain is 


1 {6»}* = (13 + 2 A)/29 ~ 0.89495, (29a) 


r _-,* V / 5[-127- 841og(x a )] 

2{A --4- 

in d = 2, and 


3 


M* - 



148 


(28a) 

and the expected convergence to this result in higher d 
as x a —> 1 is depicted in Fig. [7}). In the corresponding 
(28b) limit x a —> 0 the optimal gain in d = 2 behaves as 
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2032 


2{ey * 


v '5[-127- 841og(a: 0 )] 


- 293/5 + log(x 0 ) 


1344 


\/5[—127 — 841og(x a )] 


84 


63[3 + 41og(x a )] 

I- 


(29b) 


ultimately converging to 1/3 logarithmically slowly. Con¬ 
versely, in d = 3 the gain converges to 1/4 much more 
rapidly, 

sW* - \ “ (H + log ( ;E a)) x a . (29c) 

We therefore find that for vanishingly small targets an 
optimal heterogeneous search in a random system con¬ 
figuration with uniformly distributed starting position 
is remarkably 3 and 4 times faster, respectively in two 
and three dimensions. We should stress here that this 
gain is achieved under the constraint of a conserved D , 
which means that we do not introduce any additional re¬ 
sources. This is in striking contrast to the optimization 
of active-passive intermittent strategies, where the bal¬ 
listic excursions based on active motion are by definition 
more expensive [1, '§1 l5ll - l57j . If we were to relax the con¬ 
straint on the conserved D , the gain could, in principle, 
become arbitrarily large. 

VIII. DISCUSSION AND CONCLUDING 
REMARKS 

We analyzed the kinetics of Brownian search in 
quenched heterogeneous media. Analyzing a minimal 
model system, which captures all the essential physical 
aspects of the problem, we obtained exact analytical re¬ 
sults for the MFPT of a particle to find the target in var¬ 
ious settings. We showed that the MFPT for Brownian 
search, both homogeneous and spatially heterogeneous, 
is dominated by direct trajectories. Under the constraint 
of conserved average dynamics, we proved the existence 
of an optimal heterogeneity, which minimizes the MFPT. 

We demonstrated and explained how a blind diffusive 
searcher in a spatially heterogeneous environment can 
significantly outperform the homogeneous random walk 
when the motion is faster near the target. This gain, 
which depends on the size of the target, is significant and 
persists upon averaging over the starting position, inter¬ 
face position, or even both. The enhancement is hence 
very robust. In contrast to conventional search strate¬ 
gies (intermittent or Levy-stable motion), which have 
the highest gain in lower dimensions, the heterogeneous 
search performs best in higher dimensions. Because of 
the fact that the MFPT is dominated by direct trajecto¬ 
ries and the heterogeneity does not affect the compact¬ 
ness of exploring the surrounding space—but instead acts 
by enhancing or retarding the local dynamics—it per¬ 
forms better for non-compact exploration. 


According to recent single particle tracking experi¬ 
ments in living cells the diffusivity of smaller proteins 
is faster close to the nucleus and slower in the cell pe¬ 
riphery [GPL l6lj| . In addition, even in the presence of 
quenched spatially disordered heterogeneity, which is of¬ 
ten observed in particle tracking experiments inside cells 
©HIlnHzl, the target search kinetics can be enhanced 
as well, even if the process starts from a spatially uniform 
initial distribution of the searching molecules. Inside cells 
signaling proteins are found at extremely low concen¬ 
trations, down to, for instance, a dozen of A-repressor 
molecules searching for a single target in E. coli cell, 
whose volume is about 1 n m 3 . The search kinetics is thus 
central and rate-limiting for signaling dynamics. Hence, 
heterogeneous search processes are important and rele¬ 
vant phenomena at the few-molecule level. 

At this point a few additional remarks are in order. 
The additivity in Eq. © only holds for the MFPT and 
not for higher moments. Based on our present arguments 
the additivity principle in Eq. © should hold for an ar¬ 
bitrary number of segments but a formal proof is part of 
our current investigation. If this is indeed the case, an op¬ 
timal heterogeneity function can be formally constructed 
for any situation by first taking the appropriate limits 
of small segments and afterward optimized using varia¬ 
tional methods. The underlying physical principle how¬ 
ever, will remain unchanged. This will allow the study 
of more realistic heterogeneity profiles such as those ob¬ 
served experimentally in living cells ((ini [cl . ItTI ITht . 

In the present context the optimization of diffusion het¬ 
erogeneity cannot be perc eived as a search strategy in the 
traditional sense [51 |8I. lL8l ^ [2ll . |23( | . A searcher would some¬ 
how have to know the position of the target in order to 
optimize his motion pattern according to a optimal het¬ 
erogeneity. Conversely, in a cell the position of the target 
is given and so are the cytoplasm and nucleus properties 
giving rise to a spatially varying diffusivity. However, it 
could be part of an evolutionary optimization to improve 
the search efficiency of biomolecules in cell regulatory 
processes. 

The present results can be extended and generalized in 
numerous ways, the immediate extension being the study 
of the full distribution of first passage times. Further¬ 
more, within the context of a traditional search strat¬ 
egy our findings cannot be directly used to asses the 
efficiency of a heterogeneous search with a general off- 
center target position as in jo(| . The optimization of 
the target position-averaged MFPT with respect to the 
ratio ip of diffusivities is part of an ongoing investiga¬ 
tion. Given the present results we can speculate, how- 
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FIG. 7: a) Ratio {0} = {T}/T as a function of p = D 1 /D 2 for different dimensions and target size x a = 0.1 (full) and 
x a = 0.01 (dashed), b) Optimal heterogeneity as a function of the target size, c) {(9} for the optimal strategy as a function of 
the target size. 


ever, that an enhancement might indeed be possible in 
terms of an optimal heterogeneity at least for target po¬ 
sitions not too close to the external boundary. One can 
therefore imagine that heterogeneous search strategies 
are also beneficial for computer search or stochastic mini¬ 
mization algorithms. The ideas can be generalized to dif¬ 
fusion processes in more complex disordered systems 0- 
[76|| and even anomalous diffusion processes of continuous 
time random walk type [l3. Namely, in a finite system 
the heavy-tailed waiting time density between individual 
jumps in a subdiffusive continuous time random walk is 
expected to be exponentially tempered, exhibiting subd¬ 
iffusion over a transient but long time scale, which would 
ultimately terminate with a normal diffusion regime. In 
such a system the MFPT to the target will be finite and 
it would be interesting to investigate whether and how 
the existence and properties of an optimal heterogeneity 
change in a heterogeneous system with transiently sub- 
diffusive dynamics. The most obvious extension of the 
present results, however, goes in to the direction of het¬ 
erogeneous intermittent search. Clearly, a combination 
of both could lead to a highly superior search dynamics 
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Appendix A: Computer simulation results 

Computer simulations were performed according to the 
scheme of transition probabilities outlined in Sec. II. 
10 5 trajectories were simulated for each set of x a , Xi , 
and/or xq. The simulations correspond to a discrete ran¬ 


dom walk in between spherical shells with an absorbing 
boundary at x = x a and a reflecting one at x = 1 record¬ 
ing the number of steps within each region nip as well 
as the total number of steps n to obtain the process time 
according to 

AR? A Z? 2 

t = 7^-( n i -n 2 +n) + ^jj-(n 2 -n 1 +n). (Al) 

Note that n also counts the number of steps in the in¬ 
terfacial shells. The results for various cases are plotted 
below. We find a remarkably good agreement between 
our analytical and the simulation results. The MFPT as 
a function of the interface position is shown in Fig. [5] 
Intuitively, the MFPT depends strongly on the starting 
and interface positions and can exhibit none, one, or two 
local minima as a function of Xi . Hence, it is possible for 
a given ip that two distinct interface positions lead to the 
same MFPT. It should be noted that an equal ip does not 
correspond to equal D\ and D 2 , only their ratio is fixed. 
A degeneracy of the MFPT depending on the interface 
position is thus not surprising. The simulation results 
for the Global MFPT are shown in Fig. [9] and compared 
to the analytical predictions of Section V. In the case of 
the Global MFPT as well, zero, one or two minima are 
observed depending on ip. Here we observe significant 
differences in the features of the Global MFTP at equal 
ip in different dimension. This is due to the fact that 
the statistical weight of various starting positions is dif¬ 
ferent in different dimensions and more distant starting 
configurations have a higher weight in higher dimensions. 
The simulation results for the interface position averaged 
MFPT in a random system are shown in Fig. llOl and com¬ 
pared to the analytical predictions of Section VI. In the 
case of MFPT with a random interface position, the dis¬ 
order averaged MFPT (Fig. [Toll is only weakly dependent 
on xo as long as ip is not too large. Again, very distinct 
behavior is found for different dimensions at equal p. In 
contrast to the Global MFPT this originates solely from 
the differences in the exploration of space. The results for 
the interface position averaged Global MFPT in a ran- 
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FIG. 8: Ratio 9 = T(a;o)/T°(a;o) as function of x\ for xo = 0.6 and x a = 0.1 in different dimensions. The colors depict results 
for various heterogeneities: ip =0.6 (blue), 2 (red), 10 (green), and 150 (orange). The symbols are results of simulations and 
the full lines correspond to the analytical results in Eqs. m- 





FIG. 9: Ratio 9 = T/T° as function of xi for x a = 0.1 in different dimensions. The colors depict results for various het¬ 
erogeneities: ip =0.5 (blue), 2 (red), 10 (green), and 120 (orange). The symbols are results of simulations and the full lines 
correspond to the analytical results in Eqs. (1251) . 
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FIG. 10: Ratio {#} = {T}/T(a;o) as function of xo for x a = 0.1 in different dimensions. The colors depict results for various 
heterogeneities: ip =0.5 (blue), 2 (red), 10 (green), and 120 (orange). The symbols are results of simulations and the full lines 
correspond to the analytical results in Eqs. J20]>. 
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FIG. 11: a) Ratio {0} = {T}/T as a function of ip for different dimensions and target sizes. In dimensions 2 (b) and 3 (c) the 
target sizes correspond to 0.05 (blue), 0.1 (red) and 0.2 (green), respectively. 


dom system are shown in Fig. 1111 Also here we find an 
excellent agreement between theory and simulation. We 
see that the dependence on the target size becomes more 


prominent in higher dimensions because of the different 
statistical weight of starting positions at a given distance 
from x a and and the differences in sampling space. 
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